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Abstract. We derive computable expressions of structured backward errors of approximate 
eigenelements of ^-palindromic and *-anti-palindromic matrix polynomials. We also character- 
ize minimal structured perturbations such that approximate eigenelements are exact eigenele- 
ments of the perturbed polynomials. We detect structure preserving linearizations which have 
almost no adverse effect on the structured backward errors of approximate eigenelements of 
the ^-palindromic and *-anti-palindromic polynomials. 
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1 Introduction 

A triple {X,x,y) G C x C" x C"{x 7^ 0,j/ 7^ 0) is called an eigentriple of a polynomial 
P e P™(C"^") if 

P(A)x = and y^P(A) = 0, (1) 

where Pm(C"^") denotes the space of matrix polynomials of the form P(z) — X^jLo ^''^j'^j ^ 
C"^" and is the conjugate transpose of y. The nonzero vectors x and y are called the 
right and left eigenvectors of P corresponding to the eigenvalue A, respectively. Given a 
polynomial P and a pair (A, x) of P, the backward perturbation analysis deals with finding 
minimal perturbation AP e P„j(C"^") of P so that (A, x) becomes an eigenpair of P + AP. 
If the coefficients of the given polynomial have certain distinctive structure sometimes it is 
necessary to find a minimal perturbation having the same structure as the original polynomial 
to preserve some properties (for example, eigensymmetry). 

In this paper we restrict our attention to regular matrix polynomials. We undertake 
a detailed backward perturbation analysis of ^-palindromic and *-anti-palindromic matrix 
polynomials which we define in section O These polynomials arise in many applications such 
as in the study of rail traffic noise caused by high speed trains [25l [22l [26l [TS] . Lately there 
has been a lot of interest generated into the development of structured preserving algorithms 
and the perturbation theory of palindromic polynomial eigenvalue problems [51 [SI [IHl [101 

We denote the set of ^-palindromic or *-anti-palindromic matrix polynomials by § C 
Pm(C"^"). We choose an appropriate norm |||-|||a/ on P,„(C"^"). Given a polynomial P G § 
and (A, a;) S C x C" with x^x = 1, we determine the structured backward error 77|,j(A, x, P) 
of (A, x) as an approximate right eigenpair of P G S and construct a polynomial AP g S such 
that III AP|||a/ ~ 77|/(A, x, P) and P(A)x + AP(A)a; = 0. Moreover, we show that AP is unique 
for the Frobenius norm on C"^" but there are infinitely many such AP for the spectral norm 
on C"^". Further, for the spectral norm, we show how to construct all such AP. A similar 
analysis undertaken in [3 for certain other structures. 
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We mention that structured backward error multiplied with the structured condition num- 
ber provides an approximate upper bound on the errors in the computed eigenelements. A 
detailed sensitivity analysis including explicit expression of structured condition number of 
eigenvalues of a variety of structured matrix polynomials including palindromic matrix poly- 
nomials has been investigated in [4j. Thus structured backward errors derived in this paper 
will play an important role in the accuracy assessment of eigenelements of a *-palindromic/*- 
anti-palindromic matrix polynomial computed by structure preserving algorithms. 

Due to the lack of a genuine polynomial eigensolver, the common practice to solve a 
polynomial eigenvalue problem of degree m is to solve an equivalent generalized eigenvalue 
problem of larger size. To be specific, an n x n polynomial P of degree m is converted into an 
equivalent linear polynomial = \X + Y, X € c™"^™",y G crnnxmn ^ numerically 
backward stable algorithm is employed to compute the eigenelements of L. It is shown in 
[21] that a polynomial P G Pm(C"^") can have infinitely many linearizations. In fact 
these linearizations form a vector space of dimension m(m — + m. Analyzing backward 
error of approximate eigenpair and condition number of eigenvalues of a matrix polynomial 
Higham et al. [15l [16] have determined potential linearizations of a polynomial. 

It is well known that *-palindromic/*-anti-palindromic matrix polynomials have certain 
eigensymmetry in the spectrum, and in the eigentriple as well 22, 24, 8J. Therefore to solve a 
palindromic polynomial eigenvalue problem it is very important to preserve those structures 
in the computed eigenelements. The structured linearizations which preserve the eigensymme- 
try of *-palindromic/*-anti-palindromic matrix polynomials have been constructed in [22[|24j . 
Therefore computationally it is highly desirable to identify potential structured linearizations 
which are well-conditioned. By analyzing the structured condition number, a recipe of poten- 
tial structured linearizations of a given P G § has been produced in [Ij . 

With a view to analyzing accuracy of computed eigenelements of *-palindromic/*-anti- 
palindromic matrix polynomials we follow a similar procedure as developed in [3j for a variety 
of structured polynomials including symmetric, skew-symmetric, even and odd. Indeed, we 
consider structured backward errors ??m(A, x, P) of approximate eigenelements (A, z) of P G S 
and structured backward errors ??m(A, Am-i ® x, L; v) of approximate eigenpair (A, Am-i <8) a;) 
of widely varying structured linearzations L of P, where A^-i := [A"*""'^, ■ • . , A, 1]"^ and v is 
called the right ansatz vector, see [22]. Further, we identify potential structured linearizations 
LofPfor which 77|^(A, A„i_i(g)a;, L; v) < ar^f^ (A, x, P), for some scalar a > 0. Thus we identify 
structure preserving linearizations which have almost no adverse effect on the structured back- 
ward errors of approximate eigenelements of *-palindromic/*-anti-palindromic polynomials. 
We notice that the potential structured linearizations of T-palindromic matrix polynomials 
agree with those potential structured linearizations proposed in [1| for T-palindromic matrix 
polynomials. 

The rest of the paper is organized as follows. In section [21 we review palindromic matrix 
polynomials and their spectral symmetries. In section [3] we derive structured backward 
errors of approximate eigenpairs of palindromic matrix polynomials. In section [4] we analyze 
structured linearizations of palindromic matrix polynomials and identify potential structured 
linearizations. 

2 Eigensymmetry of palindromic matrix polynomials 

A matrix polynomial P{z) — X^jlo -^"'^i ^ IPm(C"^") is called ^-palindromic or *-anti- 
palindromic if 

P*(z) = z™P(l/z) or P*(z) = -z'"P(l/z) Vz G C \ {0} (2) 

respectively, where P*{z) — J2]Lo "^j ^^'^ * ^ ^ote that denotes the transpose 

of a matrix A and the conjugate transpose of a matrix A is denoted by A^ . We denote the 
set of ^-palindromic matrix polynomials by Sp and the set of >i=-anti-palindromic polynomials 
by Sap- Unless otherwise stated we write S for both §p and S^p. Due to the structure of the 
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coefRcients, the spectrum of a *-palindromic/*-anti-palindromic polynomial inherits a spectral 
symmetry. In fact if A is an eigenvalue of P G S then 1/A* is also an eigenvalue of P. This 
eigenvalue pairing (A, 1/A*) is known as the symplectic eigensymmetry. Table [1] gives the 
eigensymmetry and structure of eigentriples of ^-palindromic and *-anti-palindromic matrix 
polynomials. 



§ 


eigenvalue pairing 


eigentriple 


T-palindromic / T-antipalindromic 


(A, 1/A) 


(A,2;,y), (1/A,y,x) 


iJ-palindromic/i?-anti-palindromic 


(A, 1/A) 


{>^,x,y), {l/X,y,x) 



Table 1: Eigensymmetry of *-palindromic/>i=-anti-palindromic polynomials. 



The results in Table [T] follow from [21 Theorem 2.1] by extending the arguments of matrix 
pencils to matrix polynomials. Note that if A = 0, that is, if is an eigenvalue of P G § then 
oo is an eigenvalue of P as well. In this paper we consider only finite eigenvalues, although 
an infinity eigenvalue can be analyzed by considering the reversal of the polynomial or by 
considering the homogeneous polynomial, see O [221 El] • 

We now show that given {X,x) G C x C" with ||x||2 = 1 and P e S, there always exists a 
polynomial AP G § such that (P(A) + AP(A))a; — 0, that is, (A, x) is an eigenpair of P + AP. 
For X G C" with ||a;||2 = 1, we define the projection .= I — xx^ . Throughout the paper, 
we follow the convention that AP G P„(C"^") is of the form AP(z) = X^Jlo z^'^Aj. 

Theorem 2.1 Let § G {§p, §ap}- Let P E E> be given by P(z) = X]j=o ^''Aj. Suppose (A, x) G 
C X C" with \\x\\2 = 1. Set r -P(A)a; and A„, [1, A, ... , A™]^. Define 

~{x^Ajx)xx" + [(XyP*rx" + e(A)'"-^' xr'^P^] , if * = T 

-{x"Ajx)xx" + [(xyP^rx" + eA™~^a;r^P^] , if * = H 

e(AAj)*, j = : (m - l)/2 if m is odd, 
e(AAj)*, i — Q : m/2 if m is even, 

where e = 1 if S = Sp, and e = -I if S = §ap- Then P{X)x + AP{X)x = and AP G §. 
Proof: The proof is computational and is easy to check. ■ 




3 Structured backward error of approximate eigenpair 

In this section we derive structured backward error of an approximate eigenpair (A, x) of a 
polynomial P G §. The backward error of an approximate eigenpair (A, x) is defined as the 
smallest perturbation by norm, AP of P such that (A, x) is an eigenpair of P + AP. Given 
^(z) = Sjlo Z'' Aj we define norm in the following manner: 

m 

II|P|||m:-(EII^.IIm)'^'' Mg{F,2} (3) 

where M = is the Frobenius norm and A/ = 2 is the spectral norm. For a variety of norms 
on P,„(C"^") see [5]. 

Recall that a matrix polynomial P is called regular if dct(P(z)) ^ for some z G C. 
Treating (A, x) as an approximate eigenpair of a regular polynomial P we define the backward 
error of (A, x) by 

i]MiX,x,P) -.^ min {|||AP|||a/ :P(A)a; + AP(A)a; = 0}, Mg {F,2}. 
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Setting r := — P(A)a;, we have the exphcit formula [3] 

r]M{X,x,P) = ||rH2/||x||2||A™||2. 



(4) 



An exphcit formula of backward error is obtained by Tisseur |27| for a different class of norms 
on P„(C"^"). See also [8]. 

Next assume that P G S is a regular polynomial. Then we define the structured backward 
error of an approximate eigenpair (A, x) by 



min ||||AP|||m : P(A)a; + AP(A)a; = 0), M e {F,2}. 



(5) 



By Theorem 12. II it is easy to see that r]f,j{X,x,P) < oo and r]M{X,x,P) < r]f^{X,x,P). 

To derive 77f^(A,x,P) corresponding to M = 2 we use Davis-Kahan- Weinberger norm- 
preserving dilation theorem (DKW in short) which we state below. 

Theorem 3.1 (Davis-Kahan- Weinberger, jllj ) Let A,B,C and D are matrices of ap- 

^ = /i and II [a C] II2 = /i- Then there exists D 



propriate sizes. Let A, B, C satisfy 



B 



such that 



A C 
B D 



fi. All D which have this property are exactly those of the form 



D = -KA"L + - KK"f/^Z{I - L"L)^^^, 



(6) 



where := {^"^l—A^ A) ^/^B^ , L :— {fj,^l—AA^) ^/^C and Z is an arbitrary contraction, 
that is, \\Z\\2 < 1. 

For a more general version of DKW Theorem see [TT] . We use DKW Theorem in the subse- 
quent development by setting Z = in ([S]) to avoid cumbersome calculations. 

Let A,„ :— [1, A, . . . A™]^, A € C. To determine structured backward error in a conve- 
nient manner we use the projection operators which were introduced in [4] to determine 
the structured condition number of eigenvalues of >i=-palindromic/*-anti-palindromic matrix 
polynomials. lis, s G {+, — } is defined by 



n±(A™) 




ti/2+1 



a/2 



1/2 n T 



V2 ' 2 

_)^(m+l)/2_|__^(m-l)/2i T 



if m is even, 
if m is odd. 



(7) 



Now we state some basic properties of n_|_ and n_ that will be used in the subsequent 
development. It is straightforward to check that the following relations hold. 



|n+(A„,)||2- ||n_(A„ 



^(m-2)/2 2re((A)J A™-J) + |A™/2|2, if m is even 
2re((ApA'"-J), if m is odd. 



|n+(A„)||2-f ||n_(A„)||2 = ||A™| 



2|in+(A™) + n_(A„)| 



2||n+(A„)-n_(A™)| 



^^(™-l)/2|^,„_,|2^ 



^(m-l)/2 



if m is even 
if m is odd. 

if m is even 



|A^ P, if m is odd. 



3.1 T-palindromic and T-anti-palindromic matrix polynomials 

Now we derive structured backward error of approximate eigenpair of T-palindromic and 
T-anti-palindromic matrix polynomials. Recall that a polynomial P{z) = X^Jlo -^j 
palindromic if = Am-j, and T-anti-palindromic if — —Am-j- The set of T-palindromic 
and T-anti-palindromic matrix polynomials is denoted by Sp and Sap respectively. 
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Theorem 3.2 Let § e {S^jSap}. Let P £ § be given by P{z) — X]j=0'^"'^J- i^^^) "with 
\\x\\2 = 1 be an approximate eigenpair of P. Set r := — P(A)a;. Then we have 

r;|,(A,x,P) = iaUX)\\r\\l+bU\)\x^r\Y', M G {F, 2} 
where a|^(A) and fo|,/(A) are given by 



odd 



A 



A^-1 



A^l 



A = -1 



A 1 



A = ±1 



A 7^ ±1 



A^l 



yap 



a|,(A), Af = 



5|,(A), M = F 



l|n.(A„)|| 



4||n,(A„)||^-3e|A'"/^|^ 
(2||n,(A„)||i-e|A™/2|2)^ 



A 



<f(A), Af 



4||n+(A„,.)+n_(A„)||^ 
4||n+(A„o-n (A,„)||^ 



odd 



|A|>1 
|A|<1 



||A„||^-4||n+(A„)+n_(A„)||^||n,(A„)||; 

l|A„||^||n,(A„)||^ 
||A„irj-4||n+(A„)-n_(A„)||^||n,(A„)|[; 



l|A„ 



l|n.(A„.)|| 



A = ±1 







A = 1 



|A|>1 



4||n+(A„)+n_(A„ 



-|A' 



|A|<1 



4||n+(A„0-n_(A„)||^-|A"'7^ 



r 4||n,(A„)||^-3e|A'"/^|^^ 
■ 

4||n+(A,„)+n-(A„)||g-|A'"''^|^ l 



r 4||n,(A,„)||^-3e|A"'/^|^^ 
■ 

4||n+(A„)-n-(A„0||^-|A'"/^|^ i 

IIA^II^ \ 



where s = +,e — I ifS — §p, and s — —,e — ifS — Sap- 
Proof: First suppose that m is even. Then by Theorem 12. 1[ note that there aways exists a 

^ f 



polynomial AP G § which satisfies AP(A)a; + P(A)x = 0. Consequently 77|^(A,x,P) < oo. Let 



Q — [x, Qi] be a unitary matrix, where x G C" is given and Qi G C"^^" is an isometry 
so that QiX = 0. 



Let § = Sp. Define 



AAj := Q 



a(„i/2)(m/2) a„/2 



^/2 



./2, 



j = : (m - 2)/2 (8) 



and AAm-j = (AAj)^. Now AP(A)a; + P(A)x = yields AP(A)x = -P(A)a; = r(say). 
Therefore by ([8]) we have 

Y^(m-2)/2 , J, Y^(m-2)/2 \m-j„ , \m/2„ , 

By Lemma lA.2l the minimum norm solution ofJ^l^Ho^''^^ A-'5j+^j™Q"^-'^^ A™^-'aj+A™/^am/2 
Qf'r is given by 





■ T " 




a; r 








_Qfr 



|A™|li 



(A) 



lA^lli 



)f r, J = : (m - 2)/2, a„ 



./2 



(A) 



m/2 
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Further, by Lemma lA. 11 Hlthe minimum norm solution of X^jLo ^''^jj ~ given by 



if A = ±1, 



S^y+l^r^x^r ifA^il, 



a(m/2)(m/2) 



2||n+(A„)ll^-|A"/2|2 

where j — : {m — 2)/2. Thus we have 



jjx r 



2||n+(A,„)||i-|A'"/2|2 



if A = ±1, 
if A 7^ ±1, 



whenever A = ±1, and 



T 

A^ nT„ 



iff 



j = 0:m/2,AA„_, = (AAj)' 



(9) 



(A)"+(A) - - T„ l-^; ' ' (nTr\ 

n+(A„)||^-|A"/2|22; iPOlf ^'^1^'' 



(A)" 



A^. 



,/2 



(AF 
JA)™/^ 



Q^, j = 0: (m-2)/2, (10) 



2||n+(A„)||i-|A"'/2p 



V2 



(A)' 



V2 



,/2 



whenever A 7^ ±1. Setting Xj 0, j — : m, and using the fact that 
we obtain 



(11) 



(12) 



if A = ±1 



Mli + ( 



4||n+(A„ 



-3|A"'/2|2 



(2||n+(A„)ili-|AW2|2) IIAmll^ 



|a;'^r|2 if A 7^ ±1. 



Now we derive 1^2'' {X, x,P), by using DKW Theorem. If A = ±1, by ^ and Theorem 13.1 



we have /iaa, = II A A 



J 112 - Wr, 



r, j = : m, given by 



A^-^^r Qrr(gf r)^ 
' l|A™||l(||r||2-|a:^r|2)'^™/2 



A^/^x^r g?^r(g^r)^ 
l|A™||i (||r||i-|a;^rP) 



(13) 



where j = : (m - 2)/2, and X,„_j = Xf. If A 7^ ±1, by HU]), (HB) and Theorem O we 
have 



\AA 



ji|2 



|AJ+A-"-J|2 |x^r|2 , |A-"--^|2 (||^||2_|^r^|2) . , , , , , 

,„ o ,.___,„,„\^ + ' " i-^i > 1 



2||n+(A„ 



-|A™/2|2) 



/ lAJ+A— J|2 |x^r|2 |AJ|2 i\\r\\l-\;^7\^ , , , ^ . 

V (2||n+(A„oili-|AW2|2)^ ' II- 



IJ'AA„ 



given by 



= IIAA 



m/2ll2 



|A"/2|2 |2.T^|2 



|A"/2|2 (||r-||2 _ \xTr\^) 



X„ 



I (2||n+(A„)||2-|AW2|2)^ 

(|A^|^(A)'"-^+|A'"-^|^(A)^" QTrjQlrr 

(2||n+(A„)||2-|A'"/2|2) |A™-J|2 (||^j|2_|^T^|2) 

(|AJ|"(A)'"-^ + |A^'"-J|^(A)^ ^ Qfr(Qfr)^ 

(2||n+(A„)||i-|A"/2|2) |AJ|2 (||r||2_|^T^|2)' 

(A)"/2 ^ QfriQjrf 



|A™||! 



, if|A|>l, 
if |A| <1, 



72 



(2||n+(A„0||2-|AW2|2) (||r||2-|x^r|2)' 



(14) 
(15) 
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where j — : {m — 2)/2. Consequently we have 1^2^ 2;, P) — if A = ±1, and 



4||n+(A„.)+n_(A„)||^-|A"'/^P || ||2 , r 4||n+(A„)||^-3|A"'/^|^ 
« "'"^ + ^(2l|n,(A„)||i-|AW.|.)^ 

_ 4||n,(A„)+yA^)llMA-/Y ] |^T,|2 if |;,| > 

4||n+(A„.)~n_(A,„)||g-|A'"/^|^ | '||2 , r 4||n+(A„)||^-3|A-/"|" 

ll'^ll^ + L(2||n+(A„)|||_|w=p)^ 
- "'""^^'"^""a11||^"^""'"'^'1 N'^H' if |A| < 1, a ^ ±1. 



Note that if |a;"^r| — ||r||2, then ||(5^r||2 = 0. In such a case, considering Xj = we obtain 
the desired results. 



Next let § = Sap. Define 



AAj := Q 



, AA„/2 := Q 



-ai 



,/2 



O-m/2 Xm/2 



j=Q:{m- 2)/2 (16) 



and AAjn-j = — {AAj)'^ and Q = [x, Qi] is a unitary matrix. Therefore we have 



Y"(m-2)/2 , , _ ^(m-2)/2 ^m-1 

z^j=o ^ "jj Z^i=o ^ "jj 

Y^(m-2)/2 ,„i/2 Y^(m-2)/2 ,,„_j 



X r 



Note that ajj = whenever j — m/2, since (^m/2)^ = ~^m/2- By Lemma lA. 21 the minimum 
norm solution of X^j^o"^^^^ A^&j + a,„/2A™/^ — X^j^cT^^^^ X"^^^aj = Qjr is given by 



|A 



11-' = -- 



m 1 1 2 



P 1 1 2 



J'm/2 



A IP 

M i-m 1 1 2 



,j = 0: (m-2)/2. 



Also note that x^r = if A = 1. Therefore by Lemma I A. 11 [SJ the minimum norm solution of 
ElL7"^' A%,, - E\Zo'^^' A™-^a,, = x^r is given by 



[ 2||n_(A„)|| 

where j — {m ~ 2)/2. Therefore by p6)) we have 



if A = 1 
ra;'^rifA 7^ 1 



Q 
Q 







A,^ 

7 A-i 







2||n_(A„)||^2; r 



IIA^II^ 



A, 



if A = 1 



Q^if A ^ 1 



(17) 



A^„/2 = Q 







(A)'"/''(Qrr)^ ' 

A„ 



Q'^ifA^ 1, 



(18) 



where j — : (m — 2)/2 and AA„i_j = (AAj)"^. Now setting Xj — 0, we obtain the desired 

result for M = F. Further, by employing DKW Theorem l3.1l and following a similar arguments 

s 

as that in the case of § = Sp, we obtain 772°" (A, a:, P). 



Next consider m be odd. Let S = Sp. Define 



AA, := Q^AAQ = 



and (AAj)'^ = AA„,-j,j = : (m - l)/2. 
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where Q = [x, Qi],Qi G C"^^" is a unitary matrix defined as above. Consequently we 
have 

^(™-l)/2 ^ 

Note that x'^r = if A = — 1. Then by Lemma IA.2I and Lemma lA.li [3] the minimum norm 





■ T ' 




X r 








_Qlr 



solutions of Yl^jZo ^^^"^ ^^^j + Y^jZ^ '^"^ "'^^j Qi^ ^nd X^Jlo "^^^n — ^'^'^ ^"^^ given by 



^(m-l)/2 



A l|2 



A l|2 



(5j_ r, fljj — 



0_ 



if A = -1, 



^pj^(^-^x^rifA^-L 



Thus we obtain 



Q 



X 



2||n+(A,„)||^ ^ ^TC^^i 



if A 



||A,„||i^l' 

Now setting Xj ~ 0, j = : m and by we obtain 



Q^, if A 7^ -L 



/m+l ' 



l|A™|| 



if A = -1 

||A„||^||nH.(A)||^ F ^1 ,ifA^-l 



(19) 



Moreover by DKW Theorem, (|19p and following a similar techenique used for even m, we 
obtain 



4||n+(A„)+n_(A„)|| 



||A„|||-4||n+(A„)+n_(A,„)||H|n+(A„)||i|^-r 



l|n+(A,„)|| 



x^r\-, 
if|A| > 1 



4||n+(A„)-n_(A„)||| |l II2 

l|A™|li Il'^ll2 



||A,„|||-4||n+(A,„)-n_(A„ 



llA,„|i||in+(A„ 



|n.(A^)ili , TH2^ 



if |A| < 1. 



Hence the result follows when S = §p and m is odd. Following a similar arguments we obtain 
the desired results for § = Sap- B 

Remark 3.3 Observe from the above proof that ry|n(A,x,P) is obtained by the only choice 
Xj = 0. For ri2{X, x,P), by DKW Theorem, the choice of Xj is infinite. Therefore the minimal 
structured perturbation is unique for Frobenious norm and in contrast we have infinitely many 
minimal structured perturbations for spectral norm. 

Let P g §. Treating (A, x) € C x C" with ||a;||2 = 1 as an approximate eigenpair of P, we now 
construct a minimal structured perturbation AP by simplifying the expressions of AAj given 
in the proof of Theorem l3.2l Let Px I — xx^ where / is the identity matrix of order n and 
7^ a; e C". Define 



(A)^' + e(A)'"-^' 



■ 2||n,(A„)|| 



ix^r)xx", G 



{Xy + e(A)'"-^' 



;{x'^r)xx" + Fj 



F, 



1 



' - 2||n.(A™)||i-a|AW2|2V- 
™-.^.T„i „ (Xr^'ix^r)xx" 



2 + -^'"/2 



(|A^f (A)"'-J>e|A'"-Jf (Ap) xTrP];rr'^Px 
\\r\\l~\x'^r\^ 



, Lj 



(A)'"-J x'^rP];rr^Px 



where s e {+,-},£ G {+l,-l},a e {0, 1} and j e {0, 1, . . . ,to}. 
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Corollary 3.4 Let S G {Sp,Sap} cind P G S. Let (A, a;) be an approximate eigenpair of P. 
Then the unique structured perturbation AP G S when M = F, and a structured perturbation 
AP G § when M = 2,ofP for which P(A)x + AP(A)a; = and |||AP|||m = r/|^(A,a;,P) are 
given by 

L m is odd: 









AA, 


Fj i/A = -1,§ = §p 


Fj i/A -1,S = Sp 


Ej + Fj if X ^ = §p or 

X^hS^Sap 


+ F, - */l^l>l 
and S G {Sp, Sap} 


II A™ III ] 
i/ A = 1,§ = Sap 


Fi + F, - 2\\ZZ)\\i^^ 1^1 ^ 1 
A 7^ -1, § = Sp or A 7^ 1, § Sap 


[rx"" — xr^ ] if X = 1,S = Sap 



2. m is even: 





IIMII = IIMIIf 


IIMIHIIMIb 


AA, 


l|A„||^ -^-^ ^ -^J 

A = ±l,S = Sp 


a"== ±1,§ ^ Sp 




«/A7^±l,S = Sp 

or X ^ 1,S = Sap 


^'J 2||n,(A„)|||-Q|A™/2p-^J 
lf\X\ > l,Se{Sp,Sap} 






^3 2||n,(Aj)||i-a|AW2|2^3 */ |A| < 1 and 
A 7^ ±1, S = Sp or A 7^ 1, S = Sap 




z/A = l,§ = §,p 


[rx" - afr^] A = 1, S = S,p 




(x'-'r) — H , 
l|A™|l2^^ +^m/2 

i/A = ±l,S = §p 


(a'-'V) — H , T 
||yY^||i2:a; + i'm/2 

if X = ±l,S = Sp 




ff„/2 */A^±l,§ = Sp 


A,„ \\i 

F^m/2 2||n,(A„)||2-Q|A"'/2|2 -^m/Z 

if X^±l,S = Sp 




— ^^[rx^ — xr'^] if 

XeC,S = Sap 


■jp^ — p[ra; — xr J ij 

AG C,S = Sap 



i/;/ieree-I^ */§ = Sp: */S = ^' a-/^ if m is odd, 

'"'^^'^^ 1-1 */S = Sap, 1 - «/S = Sap, " " \ 1 if mis even. """^ 

Proof: First consider S = Sp. Let m be even. If A = ±1 then simplifying ^ wc have 



AA, 



{x^r)xx" + Tr^-^[\'QiQlrx" + A™-^xr^QiQf ] + Qi^^Qf 



A l|2 



A l|2' 



{x'^r)xx" + ^^[A^Pjrx^ + X'^-'xr^P,] + Q^XjQ^ 



A 112 
A^ 



1A™||1' 



{x'^r)xx" +Fj+Q^XjQ{ 



7^ V n,H 



and if A 7^ ±1 then simplifying pO]) and (|TT|) we have 

(A)^ + (A)™-^' 



AA,- 



AA 



ml'2 



2|in,(A,„2lli-|A™/2p 

(A)"'/2 
2||n,(A„0||i-|A™/2|2 



{x'^r)xx" + + QiXjQf 
{x'^r)xx" + F„,/2 + QiX,„/2Qf . 
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Now setting Xj = for M = F we obtain the unique polynomial AP G §, and putting 
Xj given in (IT31)-([IS1) for M = 2 we obtain AP G § such that P(A)a; + AP(A)x = and 
III AP|||m = '7m(A, X, P). The proof is similar when m is odd, S = Sp, and S = Sap, rn is either 
even or odd. H 

Note that if y e C"^" is such that Yx = and Y^x = then Y = {I - xx"Y Z{I - xx") 
for some matrix Z. Hence from the proof of Theorem l3.2l and Corollarv l3.4l we obtain that if K 
is a T-palindromic (resp. T-anti-palindromic) polynomial such that P(A)a:: + K(A)a; = then 
K(z) = /\Y'{z) + {I — xx^Y''ii{z){I — xx^) for some T-palindromic (resp. T-anti-palindromic) 
matrix polynomial N, where AP is given in CoroUarv 13.41 



3.2 if-palindromic and i?-anti-palindromic matrix polynomials 

We now consider the set of i?-palindromic and 77-anti-palindromic polynomials denoted by 
§p and Sap, respectively. To derive the structured backward error of approximate eigenpair 
of a polynomial P G S, where S G {§p,Sap}, we proceed as follows. Let z G C. Let us define 
the maps vec : C ^ and M : C ^ M?'^'^ by 



vec(z) 

Then we have the following Lemma 
Lemma 3.5 Let z G C and S = 



re(z)' 
im(z) 



andM(z) 



re(z) — im(z) 
im(z) re(z) 



(20) 



Then the following hold. 



1 
-1 

(i) vec(z) = Svec(z). {ii) vec(ziZ2) = ^/(zi)vec(z2), zi,Z2 G C. [iii) M{z) = M{z)'^ . 



Proof: The proof is obvious. ■ 

Note that P is iJ-palindromic polynomial if and only if iP, i := \/—\ is iJ-anti-palindromic 
polynomial. Thus the map _ff-palindromic i— s- _ff-anti-palindromic is an isometric isomorphism. 
Also observe that ??|,/(A, x, P) = ?7m(A, iP) where i§ := {iF : P G §}. We denote the Moore- 
Penrose pseudoinverse of a matrix Ahy . 

Theorem 3.6 Let P E Sp be given by P{z) = J^'jLo ■ (A, a;) with \\x\\2 ^ 1 be an 
approximate eigenpair of P. Set r — P(A)a; and A,„ = [1, A, . . . , A™]"^. Then we have 



r?7(A,x,P) 



%"(A,x,P) 



where e 



V2||r||i-|r^xP< V277(A,x,P), ^f \X\ = 1 

(2Flli-^|e^„/2)-K/P) + 2 ''^''iai?^ ^/I^l^l- 
?7(A,a;,P) 



*/|A| = 1 



(2||r| 



2 '^l'^^/2) + l''l' 



(4||n+(A,„)+n_(A„)|||-e|A"'/^P)(||r|||-|3:«rP) 

*/|A| > 1 



(2plli-^|efw2)+/P) + 



(4||n+(A)-n_(A„)||^-e|AW2|2)(||^||2_|^H^|2) 

z^A| < 1. 



if m is odd, 

1 if m is even, 



^^(m-i)/2]^ vec(a;^r) 



H(rn-2)/2 ^^m/2j 



ifm is odd, 
vec{x^r) ifm is even. 



re (A^) + re (A™-^) -im (A^) + im (A""^') 
im (AJ) + im(A™-^) re (A^) - re (A™-J) 



J = 



: (m — l)/2 ifm is odd, 
: (m — 2)/2 ifm is even, 



H 



m/2 



re(A™/2)' 
im(A™/2) 



whenever m is even, and e, is the j-th column of the identity matrix. 
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Proof: First suppose that m is even. By Theorem 12.11 it is evident that there exists a 
polynomial AP G §p for which AP(A)a; + P(A)a: = 0. Let Q = [x. Qi] be a unitary matrix 
where x is given and Qi g ^"xCn-i) jg an isometry such that QiX — 0. Define 



AA, := 



^33 



H 



a(m/2)(m/2) a„/2 



■i/2 



■/2. 



Q", j=0: (m-2)/2 (21) 



and AA^_j = {AA^)" . Since AP(A)a; + P(A)a; = 0, we have 



(m-2)/2 



Em 1 2 
J=0 ^ 







x^r 









The minimum norm solution of X)7=o ^^''^ '^"'^J "t- X)7=o ^ Q?'" is given by h 



11^ ip, flj = lA^T'"- Note that for |A| = 1, we have x^r ~ {X)"^x^r. Hence the minimum 



_ (A)" 



II A F 

I'lm II2 

norm solution of X^Jlo ^ 
have 



x^r is given by aj 



(A) 



Yxx^r, j — : m. Therefore we 



AA, 



Q 



l|A™|| 



r{Q?r)" 
X, 



Q«,AA, 



72 



.1/2 



II A™ 



AA„j_j = (AAj)J, j = : (m- 2)/2, which gives, 77p'(A,x,P) 

If I A| 7^ 1, by Lemma lA.H [71 the minimum : 
by ajj = ej_^_ir, j = : m/2 where r = [Hq Hi 



t/2 



(22) 



m+1 V II 11^ I I 

If |A| 7^ 1, by Lemma lA.H [71 the minimum norm solution of X]j=o ^''^jj ~ a;^r is given 

H(,n-2)/2 ffm/2]^ vec(x^r) and 



re (AJ) + re (A""^) -im (A^) + im (A™-^)' 
im (A^) + im(A"-J) re (A^) - re (A"-^) 



,j = 0: (m-2)/2,i7„ 



,/2 



re(A"/2)- 
im(A™/2) 



Therefore we have 



A^, 



Q 



Am 11?" 



(A) 



Q^,AA„,/2 = 



{m/2) + l ||A, 
l|A„JI? 



72 



(23) 



and AA„_j = (AAj)^ , j = : (m - 2)/2. Setting = 0, j = : m/2 we obtain 



r?F(A,a:,P) 



\nr\\i 



l'^(m/2) + l''l^) 



A l|2 



Next we derive the result for spectral norm. Note that for |A| = 1, by (|22| and Theorem l3.1l 



we have haa. 



1|A„ 



and 



A™-J r^x Qf r(gf r)^ _ A'"/^ r^x Qf r(Qf r)^ 



I|A™||1 (||r||i-|r^:rP) 



(24) 



Thus we have 772'' (A, x, P) 



= . If |A| 7^ 1 by (Hg) and Theorem [HI] we have 

if |A| >1, 
if |A| < 1, 



eT^,f|2 + ^:riM^™, if|A|>i 



|A"P 



ia;^ 

TiT^ 



lA^ 



^fm/2) + l'"l^ 



J — — given by 



|A-P (lk||^-|a;«rP) 



if |A| > 1 

if |A| < 1, 



(25) 



11 



for j = Q : {m — 2)/2. This gives 



(2||r||i 



''(m/2) + 



(4|[n+(A,„)+n_(A„ 



1|A„||^ 



772''(A,X,P)= <^ 



if|A| > 1 



(2||^ 



'("1/2) + 



1?^P) + 



(4||n+(A,„)-n_(A„ 



.|A"i/2|2)(||,.||2_|^Hr|2) 



if|A| < 1. 



Note that if \x^r\ = \\r\\2, then |lQ^r||2 = 0. In such a case, considering Xj = we obtain 
the desired results. Hence we are done for even m. FoUowing similar arguments the desired 
result can be obtained whenever m is odd.H 

Remark 3.7 Observe from the above proof that rjp {X,x,P) is obtained by the only choice 

Xj = 0. For ri2'' {X,x,P), by DKW Theorem, the choice of Xj is infinite. Therefore the mini- 
mal structured perturbation is unique for Frobenious norm and in contrast we have infinitely 
many minimal structured perturbations for spectral norm. 

Let (A, a;) with ||x||2 = 1 be an approximate eigenpair of a polynomial P e Sp. Now we 
construct a minimal structured perturbation AP by simplifying AAj given in the proof of 

pi^[A'"-^xr^P.+ 



Theorem l3.6l We proceed as follows. Let Ej 

(XyP.rx"] and K = n^l-'g;,. . 



<^{r"x)xx^, F, 



i|A„|| 



Corollary 3.8 Let P G Sp and (A, x) be an approximate eigenpair of P. Then the unique 
structured perturbation AP G S when M — F, and a structured perturbation AP when ill = 2, 
of P for which P(A)a:; + AP(A)a:; = and |||AP|||j\/ = 7]j(J(A, a;, P) are given by 





IIMII = IIMIIf 


IIMII = IIMII2 




Ej+Fj tf\X\ = l 
ej^^?xx" +Fj i/|A| 7^ 1 


E^ + Fj - x^rX'^-^WKraWt^K if \X\ = 1 
e^^^rxx" + Fj - e'^^ir(Xy X'^-'\X"'-^\-^ K if \X\ > 1 
e^^^rxx" + F, ~ e^^^?{Xy X'^-3\X'\-'^k if \X\ < 1 



where j = : (m — l)/2 if m is odd, and j — : m/2 if m is even. 

Proof: Setting Xj = when M = F and putting Xj given in (l24l)-(l26l) when M = 2 the 
proof follows by simplifying (P^ .B 

It is evident from the proof of Theorem 13.61 and by Corollary 13.81 that, if K is a H- 
palindromic matrix polynomial such that P(A)a; + K(A)a; = then K(z) = AP(z) + {I — 
xx^)N{z){I—xx^) for some iJ-palindromic polynomial N, where AP is given in Corollarv l3.8l 



4 Structured backward error and palindromic lineariza- 
tions 

The classical way to solve polynomial eigenvalue problem is to convert the polynomial P into 
an equivalent linear polynomial L, called a linearization of P, and compute the eigenelements 
of L. For a polynomial P the set of linearizations form a vector space Li(P) defined by [22] 

Li(P) :={L(A) :L(A)(A„_i®/) = w®P(A)},i;eC™, (27) 

where v is called the right ansatz vector, (g) is the Kronecker product, A,„_i := [A™~^, ... A, 1]"^ 
and L is of the form L(A) = XX + Y,X e C""^"", F g ^-mnxmn arbitrary lineariza- 

ton L G Li(P) can destroy the eigensymmetry of a structured polynomial P [52]. Hence to 
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solve a structured polynomial eigenvalue problem one needs to choose a linearization which 
preserves the eigensymmetry of the polynomial. These linearizations are called structured 
linearizations. 

Mackey et al. [25j have shown that a *-palindromic/*-anti-palindromic matrix polynomial 
can have both *-palindromic and *-anti-palindromic linearizations that preserve the eigen- 
symmetry of the polynomial. Table [2] gives the structure of ansatz vectors for structured 

1 " 

for details see [2l [25] 



linearizations, where R = 



1 





Structured Linearization 


ansatz vector 


T-palindromic 


T-palindromic 


Rv = V 


T-anti-palindromic 


Rv —V 


T-anti-palindromic 


T-palindromic 


Rv = —V 


T-anti-palindromic 


Rv = V 


_ff-palindromic 


_ff-palindromic 


Rv = V 


_ff-anti-palindromic 


Rv = —V 


_ff-anti-palindromic 


i/-palindromic 


Rv = —v 


_ff-anti-palindromic 


Rv — V 



Table 2: Table of the admissible ansatz vectors for palindromic polynomials. 



Note that to solve palindromic polynomial eigenvalue problem, the prime task is to de- 
tect potential structured linearizations that behave well during computations. Analyzing 
sensitivity of eigenvalues, potential structured linearizations have been produced in j3] for 
T-polynomial/T-anti-palindromic matrix polynomials. With a view to analyze accuracy of 
approximate eigenelements, in this section, we identify the potential structured linearizations 
of a *-polynomial/*-anti-palindromic polynomial. 

We first review some basic results available in the literature. Let P be a regular polynomial. 
Let L G Li(P) be a linearization of P corresponding to the right ansatz vector v G C™. Then 
the relationship between the eigenelements of P with that of its linearization L is given in 
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• X e C" is a right eigenvector for P corresponding to an eigenvalue A € C if and only if 
Am-i (E) a; is an eigenvector for L(A) corresponding to the eigenvalue A. 

Treating (A, x) as an approximate eigenpair of P the relations 

||L(A)(A„_i®x)||2 = ||i;||2||P(A)x|l2, (28) 
|(A„_i®a;)^L(A)(A™_i®x)| = \Al_,v\\x'^P{X)x\, (29) 
|(A™_i(5§a;)^L(A)(A„,_i®:E)| - \A^_,v\\x''P{X)x\ (30) 



have been derived in [15]. In view of (p8]) -(]30 p . without loss of generality we assume that the 
right ansatz vector v is of unit norm. Note the inequality 



2m - ||A™_i||2||(A, 1)||2 - ^ ' 

given in [TBI Lemma A.l]. 

We denote the backward error of {X,Am-i ® x) by 77(A,Am_i (g) x,L;v). Comparing 
rjM{X, X, P) with TjMiX, Am_i ® x, L; v) the inequality 



2m ?7M(A,a;,P) 
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has been proved in [7, Theorem 4.1]. Now recaU that 77m(A, a;,P) < rif_j{X,x,P). Hence for 
any structured hnearization L g Li(P) of a given P g §, we have 7/m(A, Am_i (8) x,h]v) < 
ri%i{X, Am_i (g) X, L; v) where the ansatz vector v is given in Tabled Thus by [31 Lemma 4.2] 
we have 

m_M^ ,|,(A,A ,^.L;.) ^^^ (33) 
2m 77Af(A,x,P) 
In the sequel we use the inequahty 

0<^^<i,*e{r,i/} (34) 

||Am-l||2 

where v e C™ with ||w||2 = 1. 

Theorem 4.1 Let P be a T -palindromic matrix polynomial. Let §p C Li(P) and §ap C 
Li(P) be the set of T -palindromic and T-anti-palindromic linearizations ofP, respectively. 
Suppose Lp S Sp and L^p S Sap are t/ie T -palindromic linearization and T-anti-palindromic 
linearization of P corresponding to the ansatz Rv — v and Rv — —v, respectively. If {X,x) 
with \\x\\2 = 1 is an approximate right eigenpair ofP then we have 

• M = F : 

1, re(A) > : ■/l-^,/^< n%^{\^^-.®x,W,v) < ^ 
V V ^ - ri{X,x,P) 

rip''{X,x,Lap;v) 

• M = 2 



m + l ^ T72°''(A,X,Lap;t)) 



Proof: Let r := -P(A)a;. First consider M = F. By Theorem IX2l and using (l28l)-(l29l) we have 

l|A™||2yi|r||l-2^^-i^|a:^ 



r?/(A,a;,Lp;'i;) 



x/2 



?7m(A,x,P) IklblKl, A)||2||A™-i||: 

if A 7^ —1, and 



77/" (A, a;, Lap; u) 



^ l|A™||2yilii + 2^^g-i^|xT,|2 
\/2 ^ 



mi{\x,P) |lr||2||(l, A)||2||A™_i||, 



if A ^ 1. It is easy to verify that ^/l - < J\\r\\l^ 2^ ^g-^J jx^r^ < ||r||2 



if A ^ -1, and ^1 + < + 2pfS(^ |x^rP < ||r||2 if A ^ 1. Thus by 

(1211), (ini) and (03]) the desired result follows for M = F. 

Now consider M = 2. If |A| > 1 then by Theorem E^l (I1S1)-(I111) we have 



lff{\,Kn-l®X,hp-v) ^ V2|jA™||2 / , |A| 



77m(A,x,P) - ||A„_i||2||(l, A)||2 V |l + AP||A„_i||2 ||(1, A)||i 
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by ([M]) and Hence the results follow by Similarly if | A| < 1, A 7^ 1, by Theorem 
we have 



< V2 Jl 



Hence the result follows by ((33|) . ■ 

Note that |1 - A| < |1 + A| when re(A) > and |1 + A| < |1 - A| when re(A) < 0. 

Remark 4.2 Let P be a T -anti-palindromic polynomial. Then we obtain similar bounds from 
Theorem \4.1\ by interchanging the role of hp and Lap. 



Theorem 4.3 Suppose that the assumptions of Theorem \4.1\ hold. Let § be the set of T- 

palindromic polynomials. Then we have 

^ M^F: ^^"(^'^-^^^'L,;.) ^^^^^^ '^^^t"'^;-;"^ < V2./re(A)<0. 

ri%{\,x,P) ri%(\,x,P) 

V%{X,x,P) J ^ V%{\x,P) ^ ^ 

Proof: The proof is followed from the fact that riM{X, x, P) < '7|/(A, x, P) and Theorem l4.1I B 

Remark 4.4 Let § be the set of T -anti-palindromic polynomials and P G §. Then the sim- 
ilar bounds hold for T -palindromic and T -anti-palindromic linearizations but the role of T- 
palindromic linearizations and T -anti-palindromic linearizations get exchanged in Corollary 

El 



The moral of the Theorem 14.11 and Corollarv l4.3l is as follows. For a T-palindromic polyno- 
mial the bounds derived above advice to choose T-palindromic linearization when re(A) > 0, 
and choose T-anti-palindromic linearization when re(A) < 0. Observe that our choice of struc- 
tured linearizations is compatible with that given in A, by analyzing structured condition 
number. 

Now we consider i7-palindromic/i?-anti-palindromic matrix polynomials. 

Theorem 4.5 Let Sp be the set of H -palindromic matrix polynomials and P G §p. Let S C 
Li(P) be the set of H -palindromic or H -anti-palindromic linearizations of P. Suppose L G S 
is an H -palindromic linearization or H -anti-palindromic linearization of P corresponding to 
the ansatz Rv = v or Rv = —v, respectively. Lf {X,x) with \\x\\2 = 1 is an approximate right 
eigenpair ofP and |A| = 1 then we have 



1 ^ ??f (A, A,n-i (g)x;L,t;) ^ ^ r]%{\,Am-i'»x;L,v) ^ ^ 



2m r]F{X,x,P) rf/{\,x,P) 

m-\-l ^ ?7|(A,A^-i iS)x;L,v) _ 7?f(A,Am-i <»x;L,v) ^ ^ 
2m - 772 (A, a;, P) 772''(A,x,P) 

Proof: Let r := — P(A)a;. If |A| = 1, by Theorem 13.61 we have 



A II /Oil 112 |A«_it-Pb«rP 

v{\x,P) ||r||2||(l, A)||2||A™-i||2 



It is easy to verify that ||r||2 < y2||r||2 ^^^p^ — j-jp < v2||7'||2- Hence the desired result 

follows by ([55]) . For the spectral norm we obtain the desired result by noting that 77^(A, x, P) = 
rj{X,x,P) by Theorem[ 
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This shows that there is almost no adverse effect of structured linearization of a H- 
palindromic matrix polynomial on the backward errors of approximate eigenelements when 
the approximate eigenvalues are on the unit disk. On the other hand, if |A| ^ 1, a little 
calculation gives the following bounds: 



ry^(A,A ,^.;L,.) ^ ^ A + M| ^ |A| ^ 1 
'n{\x,v) y \\r\\l 



r?(A,x,P) 



< { 



where r \=r„ 



tion and r := fa 



+ re(A) im(A) 
im(A) 1 - re(A) 

"1 - re(A) -im(A) 
-im(A) 1 + re(A) 



V X 



vx 



im(A^ 

^ ^re(A^. 
im(A^ 



^P(A)x)' 
P(A)x) 



H 



for /f-palindromic lineariza- 



P(A)a:) 
P(A)x) 



for if-anti-palindromic 



linearization. 

Therefore the moral for i?-palindromic/i/-anti-palindromic matrix polynomials is as fol- 
lows. If eigenvalues are on the unit disk then it does not matter whether we choose H- 
palindromic or i7-anti-palindromic linearization. However, for eigenvalues not on the unit 
disk, it may be a good idea to solve both iJ-palindromic or i7-anti-palindromic lineariza- 
tions and pick an eigenpair (A, a;) from ff-palindromic or iJ-anti-palindromic linearization 
according as rp < rap or rap < rp. 



5 Conclusion 

We have derived computable expression of structured backward errors of approximate eigen- 
pairs of *-palindromic/*-anti-palindromic matrix polynomials. We mention that these ex- 
pressions have an important role to play in analyzing stability of structured preserving algo- 
rithms. Finally structured backward errors have been used to determine potential structured 
linearizations of a *-palindromic/*-anti-palindromic matrix polynomial. 
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A Appendix 

The following Lemma summarizes few auxiliary results required to prove the Theorems in 
section [3] 

Lemma A.l Let A e C and A„ = [1, A, A"]^ G C™+i. Let x,r e C",e G {+1,-1} 
and s G {-f,— }. Then the solution of'^"^^QajjX^ — x*r,ajj G C, * G {T, i/} that minimizes 

1. Ojj — if x*r — 0. 

n _ W * 

3. If m is odd, Ojj = f-0'(m~j){m-j)T i — ■ {m^l)/2 then ajj = ^21111^(1^'* )|p •^*^' whenever 
s = + if e = \, and s — — if e ~ —1. 

4. If m is even, Ojj = a(^ra-]){m-j), j = : m/2 then ajj = 2||nf(A'J)|||-|AW2|2 a^*'" <^nd 



a(m/2)(m/2) — 2||n4A„ ) |||- 1 A"/^ |2 ^ 



5. Ifm is even, Ujj = -a(,„_j)(„j_j) , j = : m/2 thenojj = 2||n_ (aL)II^' a"*^ a(m/2)(m/2) 
0. 



6. If m is odd, Ojj — eam-j,m-j j — ■ (™ — l)/2 and |A| 7^ 1 then ajj = eJ^iT where 
r^[Ho Hi ... H 1^,^-1)/ 2]^ yec{x*r), Hj = 

Cj is the j-th column of the identity matrix and vec is defined in I120\) 



re (A^) + ere (A^-^) -im (A^) + eim (A™"^')' 
im (AJ) + eim(A"-J) re (A^) - ere (A"-J) 
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7. If m is even, ajj = £ajn-j.m-j j — : m/2 and |A| 7^ 1 then ajj — ej^i^ where 



H, 



m/2 



■re(A'"/2) 



r^[Ho Hi ... i?(,„-i)/2 vec(x*r), 

if e = — 1 , and Hj , j = : 



"0 -1' 




■re(A™/2)" 


1 




im (A'"/2) 



if f = 1, Hra/2 = 

(to — 2)/2 is same as that given 
where Ils{Am) is defined in 

Proof: The proof of[T]and[2]are obvious. Now consider [31 Let m be odd. Then 

m (m-l)/2 



3=0 j=o 



X r. 



Hence the result follows by [2l The proof is similar for |4] and [5] To prove [6] we proceed as 
follows. Let e = 1. Apply the map vec at both sides of X^Jlo ^^'^jj ~ ^*''"- This yields 



vec(^ A*ajj) = vec(x*r) ^ ^ M(A*)vec(ajj) = vec(a;*r) 
3=0 3=0 



(35) 



where vec and M are defined in (|20p . Employing the condition ajj — am-j,m-j j — ■ (to— l)/2 
on ([35|l we have 

(m-l)/2 

^ (M(A^) + M(A™-^)E)vec(ajj) = vec(a;*r) 
3=0 



where S = 



1 

-1 



Thus we have ajj — ^J+i^ where 



[Ho Hi ... H(^„,-i)/2]^ vec{x*r), Hj 



re (A^) + re (A™-^) -im (A^) + im (A™-^)' 
im (A^) + im(A'"-^) re (A^) - re (A"-^) 



Cj is the j-th column of the identity matrix. The proof is similar for e = — 1 and[7lB 

Lemma A. 2 Let A e C, Xj E C"^^, j = : to. Then the solution ofJ2Y=Q •^"'^j = VtV ^ 
that minimizes X]j=o ll^jlli given by Xj — -ji^^y where A™ = [1, A, . . . , A™]-^. 

Proof: The proof follows by using Moore-Penrose pseudoinverse of A^.B 
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